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ABSTRACT 



Archival XMM-Newton EPIC observations reveal the 
population of X-ray sources of the bright Local Group 
spiral galaxy M 31, a low- star-formation-rate galaxy like 
the Milky Way, down to a 0.2-4.5 keV luminosity of 
4.4 X 10^^ erg s~^ . With the help of X-ray hardness ratios 
and optical and radio information different source classes 
can be distinguished. The survey detected 856 sources in 
an area of 1.24 square degrees. Sources within M 31 are 
44 supernova remnants (SNR) and candidates, 18 super- 
soft sources (SSS), 16 X-ray binaries (XRBs) and can- 
didates, as well as 37 globular cluster sources (GIC) and 
candidates, i.e. most likely low mass XRBs within the 
GIC. 567 hard sources may either be XRBs or Crab-like 
SNRs in M 3 1 or background AGN. 22 sources are new 
SNR candidates in M 31 based on X-ray selection cri- 
teria. Time variability information can be used to im- 
prove the source classification. Two GIC sources show 
type I X-ray bursts as known from Galactic neutron star 
low mass XRBs. Many of the M 31 SSS detected with 
XMM-Newton, Chandra and ROSAT, could be identified 
with optical novae. Soft X-ray light curves can be deter- 
mined in M 3 1 center observations for several novae at a 
time opening a new area of nova research. 

Key words: galaxies: individual (M 31), novae, cata- 
clysmic variables, supernova remnants. X-rays: galaxies. 
X-rays: binaries. X-rays: bursts. 



1. INTRODUCTION 



In the XMM-Ne wton survey of the Local Group Sc 
galaxy M 33 (Pie tsch et al.Ll2QQ4aL hereafter PMH2004), 
408 sources were detected in a 0.8 square degree field 
combining the counts of all EPIC instruments, which 
could be identified and classified using X-ray colors and 
time variability as well as optical and radio information. 
This proved to be an efficient way to separate super- soft 
X-ray sources (SSSs) and thermal supernova remnants 
(SNRs) in M 33 from Galactic stars in the foreground 
and "hard" sources. These hard sources may be either X- 
ray binaries (XRBs) or Crab-like SNRs in M 33 or active 



galactic nuclei (AGN) in the background of the galaxy. 
The success of this survey inspired us for a similar anal- 
ysis of all archival XMM-Newton observations of M 3 1 . 



The Andromeda galaxy M 31 is locate d at a d i stance 
similar to th e one o f M 3 3 (780 kpc, Holland, 19981 
Stanek & GarnavichL Il998l i.e. corresponds to 
3.8 pc and the flux to luminosity conversion factor 
is 7.3x10^^ cm^) and - compared to the near face- 
on view of M 33 - is seen under a higher inclina- 
tion (78°). The optical extent of the massive SA(s)b 
galaxy can be approximated by an inclination-corrected 
D25 ellipse with a large diameter of 153.3' a nd axis 
ratio of 3.09 ("de Vaucouleurs et al.", '1991'; "TullV, \E 



With its moderate Galactic foreground absorption (A^h= 
7x10^° cm-2, Star k etal.1 11992). M 31 is well suited 
to study the X-ray source population and diffuse emis- 
sion in a nearby spiral similar to the Milky Way. M 31 
was a target for many previous imaging X-ray mis- 
sions. The Einstein Observatory detected 108 individ- 
ual X-ray sources brigh ter than 5 x 10^^ erg s~^ (see e.g. 
iTrinchieri & Fabbiand [l99 1). The sources were iden- 
tified with young stellar associatio ns, globular clusters 
(GCs) and SNRs. The i^OSAT HRI (Primini etaliri99l 
detected 86 sources brighter than ~10^^ erg s~^ in the 
central area of M 31. The ROSAT PSPC covered the 
entire galaxy twice in surveys conducted one year apart 
and detected altogether 560 X-ray sources down to a 
limit of ~5xl0^^ erg s~^ and SSS were established as 
a new class of M 31 X-ray sources (Supper et al., 19971 
l200l . The flux of many of the sources varied signif- 
icantly between the Einstein and ROSAT observations. 
Deep Chandra ACIS I and HRC observations of the cen- 
tral region (covered areas of 0.08 and 0.27 square de- 
gree ) resolve d 204 and 142 X-ray sources, respectively 
(Ga rcia etal.L l2000: Kongetal, 2002b; Kaaret, 200^. 
A synoptic study of M 3 1 with the Chandra HRC cov- 
ered "most" of the disk (0.9 square degree) in 17 epochs 
using short observations, and resulted in mean fluxes 
and long -term l ight curves for the 166 objects detected 
( Willia ms et al.L ^004). In these o bservations, several 
M 31 SNRs were spatially resolved ('Kong et al.L l2002al 
2003) and bright XRBs in globular clusters and SSSs 
and quasi-soft sources (QSSs) co uld be c haracterized 
(IpTstefano et aL.200Z.2004:.Greiner et allf2004) . 
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One of XMM-Newton's most important contributions to 
galaxy science was a deep survey of the central region 
around the long axis of M 3 1 as part of the guaranteed 
time program. This survey is unique in that it has the 
greatest depth (~10^^ erg s~^) and best spatial resolu- 
tion of any existing large area M 3 1 survey. It has covered 
almost 3° (>40 kpc) along the major axis of the galaxy 
and 30' (~7 kpc) along the central portion of the minor 
axis. These deep XMM-Newton observations have al- 
lowed us for the first time to study the short-term time 
variability (~100 s and shorter, see below) and spectra of 
bright X-ray sources in a galaxy o utside the Milky Way 
and the Magellanic Clouds (e.g. 'Osborne et al.^. I200lt 
[Barnard etal, 2003, 2005; Mangano et al, 2004|). The 
observations revealed diffuse emissi on from the hot ISM 
in the centre and the northern disk iShirev etaHllool 
[Trudolvubovetal, 2001, 2004), and were used to de- 
rive source luminosity distributions iTrudolvubov et all 
|2002). 

Here, I summarize results of our group mainly based on 
the archival M 31 XMM-Newton observations, includ- 
ing X-ray images and a source catalogue for the archival 
observations of M 31 (Pietsch et al., 2005b, hereafter 
PFH2005), the det ection of type I X-ray bursts in M 31 
jPietsch & HaberlL l2005». hereafter PH2005) and on the 
detection of optical novae inM31asSSSs ( Pietsch et all 
12005^ hereafter PFF2005). 



2. XMM-NEWTON SURVEY OF M 31 



PFH2005 have created merged medium and thin filter im- 
ages for the three EPIC instruments, in five energy bands 
(0.2-0.5 keV, 0.5-1.0 keV, 1.0-2.0 keV, 2.0-4.5 keV, and 
4.5-12 keV), using only times of low background from 
the archival XMM-Newton M 31 observations which at 
the time contained four observations of the centre area of 
M 3 1 separated by half a year, two pointings in the south- 
ern, three in the northern disk, and one short observation 
in the halo, which all at least partly cover the optical D23 
ellipse. In total, the observations in the analysis cover 
an area of 1.24 square degrees (see Fig.[I} with a lim- 
iting sensitivity of 4.4x10^^ erg s~^ in the 0.2-4.5 keV 
band which is a significant improvement compared to the 
Chandra surveys. However, up to now only about 2/3 of 
the optical M 31 extent (D25 ellipse) are covered with 
a rather inhomogeneous exposure. There were still sig- 
nificant offsets between the observations that had to be 
corrected for before merging. For the centre observations 
these offsets were determined from source lists of the in- 
dividual observations using the USNO-Bl, 2MASS, and 
Chandra catalogues to define an absolute reference frame. 
This finally resulted in a residual systematic position er- 
ror of less than 0.5". The source detection procedures 
revealed 856 sources using simultaneously 5x3 images 
(5 energy bands and PN, MOSl and M0S2 camera). For 
the pointings into the disk and halo of M 3 1 this proce- 
dure was applied to the individual observations. The cen- 
tre pointings strongly overlap and therefore the images 
were merged to reach higher detection sensitivity. 




Figure 1. Combined XMM-Newton EPIC image in the 
0.2^.5 keV band smoothed with a 20" FWHM Gaussian. 
Orientation and optical D25 ellipse are indicated. 
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HRl 

Figure 2. Hardness ratios (HR) detected by XMM-Newton EPIC. Shown as dots are only sources with HR errors smaller 
than 0.2 on both HRi and HRi + 1. Foreground stars and candidate are marked as big and small stars, AGN and 
candidates as big and small crosses, SSS candidates as triangles, SNR and candidates as big and small hexagons, GlCs 
and XRBs and candidates as big and small squares. In addition, we mark positions derived from measured XMM- 
Newton EPIC spectra and models for SSSs (SI to S4) as filled triangles, low mass XRBs (LI and L2) as filled squares, 
SNRs (NI32D as NI, IE 01 02.2-72 1 9 as N2, NI57B as N3, Crab spectra as CI and C2 ) as filled hexagons, AGN (AI and 
A2) as asterisk (extracted from Fig. 5 ofPFH2005). 



Table I. Summary of identifications and classifications of 
XMM-Newton X-ray sources in the M 31 and M 33 fields 
(see PFH2005 and PMH2004). 
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SNR 
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23-2 


GIG 


27 
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XRB 
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hard 




567 




267 



Hardness ratios were calculated only for sources for 
which at least one of the two band count rates had a sig- 
nificance greater than 2a (Fig. 13. In search for identifi- 
cations, the X-ray source positions were correlated with 
sources in the SIMBAD and NED archives and within 
several catalogues. The cataloged X-ray sources are 
"identified" or "classified" based on properties in X-rays 
(hardness ratios (HR), variability, extent) and of corre- 
lated objects in other wavelength regimes. A source is 
counted as identified, if at least two criteria secure the 
identification. Otherwise, it is only counted as classified. 

Table [l] summarizes identifications and classifications 



according to the XMM-Newton catalogues of M 3 1 and 
M 33. For the SNRs in M 33 two new optical counterparts 
for soft X-ray SNR candidates from the XM M-Newton 
list are indicated fsee'Ghavami an et all 12005). Detection 
of strong time variability in follow-up analysis will cer- 
tainly move many objects from the "hard" to "XRB" clas- 
sification. Comparison to earlier X-ray surveys revealed 
transients not detected with XMM-Newton, which add to 
the number of M 31 XRBs. Up to to now, only low mass 
X-ray binaries have been identified in M 31, mostly by 
correlations with globular cluster sources. In M 33, how- 
ever, besides the ultra-luminous X-ray source close to the 
nucleus (X-8, most likely a black hole XRB) the only 
other XRB identified is the eclipsing high mass XRB X- 
7 with an orbital period of 3.45 d that has been confirmed 
by the e llipsoidal heating lig ht curve of its optical com- 
panion (jPietsch etal.t|2004b|). 

Many foreground stars, SSSs and SNRs can be classified 
or identified. The number of 44 SNRs and candidates 
more than doubles the X-ray-detected SNRs. 22 sources 
are new SNR candidates in M 3 1 based on X-ray selec- 
tion criteria. Another SNR candidate may be the first 
plerion detected outside the Galaxy and the Magellanic 
Clouds. Additional SNR candidates can be identified by 
comparing cataloged X-ray sources with o ptical narrow 
filter i mages of the Local Group survey of iMassev et alJ 
( [20011 ). Figure [3] gives two examples for two new optical 
SNR candidates proposed by X-ray source position and 
a [S II]/Ha ratio of the optical emission characteristic for 
optical SNRs. 

However, besides a few clearly identified XRBs and 
AGN, and SNR candidates from positions in other wave- 
lengths, we have no clear hardness ratio criteria to se- 
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Figure 3. SNR candidates from overlay of Local Group 
survey images (Ha above, [SII] below): [PFH2005] 224 
was classified as SSS candidate, [PFH2005] 234 was un- 
classified. 
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lect XRBs, Crab-like SNR or AGN. They are all "hard" 
sources (567 sources classified as hard in total). Only 
additional criteria like short or long term X-ray variabil- 
ity or detailed spectral modeling will reveal their na- 
ture. Such methods can be used in the M 31 center area 
with four overlapping observations in the XMM-Newton 
archive (separated by half a year) and additional observa- 
tions of the area taken 2.5 yr later (see Fig.|4|and|5l). 



Figure 4. Variability of X-ray sources within four over- 
lapping XMM-Newton observations to the M 31 centre 
performed from June 2000 to January 2002. 



DETECTION OF TYPE I X-RAY BURST 
SOURCES IN M 31 



Within the Milky Way, bright globular cluster X-ray 
sources were identified as low mass XRBs. Many of 
them show type I X-ray bursts identifying them as neu- 
tron star systems. PH2005 searched for X-ray bursts in 
XMM-Newton archival data of M 3 1 sources which were 
identified or classified as globular cluster sources in the 
PFH2005 catalogue (Fig. Two bursts were detected 
simultaneously in EPIC pn and MOS detectors and some 
more candidates in EPIC pn. The energy distribution 
of the burst photons and the intrinsic luminosity during 
the peak of the bursts indicate that at least the strongest 
events were type I radius expansion burst radiating dur- 
ing maximum at the Eddington limit of a 1 .4 Mq neutron 
star for hydrogen-poor matter (Fig. |6|). Standard type I 
bursts would show harder spectra and would not be bright 
enough to be detected by XMM-iVewtoiiEPIC. The bursts 
identify the sources as neutron star low mass XRBs in 
M 31. These type I X-ray bursts are the first detected 
outside the Milky Way and show that, with the large col- 
lecting area of XMM-Newton, X-ray bursts can be used 
to classify neutron star low mass XRBs in Local Group 
galaxies. 



M31 burst source [PFH2005] 253 
XMM-Newton, EPIC 0.2 - 7.0 keV 




4.38x10^ 4.4x10* 
Time (s) 



Figure 6. Combined XMM-Newton EPIC light curve of 
a type I X- ray burst of source JPFH2005] 253 in M31 
(Fig. 3 from\Pietsch & HaherlUcM). 
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Figure 5. Variability of X-ray sources between combined XMM-Newton M 31 center observations of June 2000 to January 
2002 (left, globular cluster sources and candidates marked, burst sources numbered) and observations 2.5 yr later ( right). 
New transients are marked. However, many bright sources from left image are missing. 



4. OPTICAL NOVAE AS MAJOR CLASS OF SSS 
INM31 



PFF2005 searched for X-ray counterparts to optical no- 
vae detected in M 31 and M 33. They combined an M 31 
optical nova catalogue from the WeCAPP survey with op- 
tical novae reported in the literature and correlated them 
with the most recent X-ray catalogues from ROSAT, 
XMM-Newton and Chandra, and - in addition - searched 
for nova correlations in archival data. They report 21 X- 
ray counterparts for novae in M 31 (mostly SSS). Their 
sample more than triples the number of known optical 
novae with super-soft phase. Most of the counterparts 
are covered in several observations which allows to con- 
strain X-ray light curves of optical novae (see Fig. Q. 
Selected brighter sources were classified by their XMM- 
Newton EPIC spectra. Six counterparts are only detected 
in Chandra HRC I (3) or ROSAT HRI (3) observations, 
i.e. X-ray detectors with no energy resolution, and there- 
fore can not be classified as super-soft. From the well- 
determined start time of the SSS state in two novae one 
can estimate the hydrogen mass ejected in the outburst 
to ~ 10~^Mq and ~ 10~^Mq, respectively. The super- 
soft X-ray phase of at least 15% of the novae starts within 
a year. At least one of the novae shows a SSS state last- 
ing 6.1 years after the optical outburst. Six of the SSSs 
turned on between 3 and 9 years after the optical discov- 
ery of the outburst and may be interpreted as recurrent 
novae. If confirmed, the detection of a delayed SSS phase 
turn-on may be used as a new method to classify novae 
as recurrent. At the moment, the new method yields a 
ratio of recurrent novae to classical novae of 0.3. Ongo- 
ing optical and X-ray monitoring of the central region of 
M 3 1 , where most of the novae are detected, should allow 




200 400 600 800 1000 
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Figure 7. Light curves for M 31 and M 33 novae that were 
detected within 1000 d after outburst. Detections of indi- 
vidual novae are connected by solid lines, and connec- 
tions toj^££er_lirmtsa^ by dashed lines ( Fig. 3 
from lPietsch et at \ \2005d) . 

us to determine the length of the plateau phase of several 
novae and, together with the nova temperature develop- 
ment, give a handle on the masses of the white dwarfs 
involved. 



5. CONCLUSIONS 



The sensitivity of XMM-Newton and Chandra obser- 
vations of M 31 combined with the wealth of multi- 
wavelength data for the galaxy allows a detailed study 
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of the point source population. Many more interest- 
ing results can be expected from further monitoring of 
M 3 1 with XMM-Newton and Chandra specifically also 
in the energy band below 0.5 keV. The first light curves 
of the SSS state of optical novae proved that these kind 
of studies can be more efficiently achieved by observing 
many candidates at the same time in one field in M 31 
than by monitoring individual novae in the Milky Way 
or the Magellanic Clouds. The results of the Chandra and 
XMM-Newton observations of M 3 1 demonstrate the im- 
portance of arcsec spatial resolution, broad energy cover- 
age, good energy resolution, and high collecting power - 
used together with deep images and catalogues at other 
wavelengths - also for future X-ray source population 
studies in nearby galaxies. 
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